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Abstract 

Using a pulse-fit method, we investigate the spectral lags between the traditional gamma-ray band 
(50-400 kcV) and the X-ray band (6-25 keV) for 8 GRBs with known redshifts (GRB 010921, GRB 
020124, GRB 020127, GRB 021211, GRB 030528, GRB 040924, GRB 041006, GRB 050408) detected 
with the WXM and FREGATE instruments aboard the HETE-2 satellite. We find several relations for 
the individual GRB pulses between the spectral lag and other observables, such as the luminosity, pulse 
duration, and peak energy -E pca k- The obtained results are consistent with those for BATSE, indicating 
that the BATSE correlations are still valid at lower energies (6—25 kcV). Furthermore, we find that the 
photon energy dependence for the spectral lags can reconcile the simple curvature effect model. We discuss 
the implication of these results from various points of view. 
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1. Introduction 

Gamma- Ray Bursts (GRBs) are the most energetic ex- 
plosions in the universe. Past studies have found that 
GRBs consist of ultra-relativistic outflows with collimated 
jets at cosmological distances. However, it is not clear how 



the central engine forms and how the electrons or pro- 
tons are accelerated in shocks and photons are radiated. 
In addition, GRBs are quite important as candidates for 
distance-indicators. Owing to their very intense bright- 
ness, GRBs can be a powerful tool to measure distances 
in the high redshift universe. 
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One of the characteristics of GRB prompt emission is 
the spectral lag, which is the time delay in the arrival of 
lower-energy emission relative to higher-energy emission. 
The previous analyses have been done using a sample of 
many BATSE GRBs between typical energy bands 25—50 
keV and 100—300 keV, using both CCF (cross correlation 
function; e.g., Norris et al. (2000) ) and peak-to-peak dif- 
ference (e.g., Hakkila et al. (2008) ). An anti-correlation 
between the spectral lag and the luminosity exists for the 
BATSE GRBs above 50 keV energies. Since we can obtain 
the intrinsic luminosity of GRBs from the lag-luminosity 
relation once we measure the spectral lag, the distance 
of the GRBs can be derived from the observed flux. But 
it is not clear whether the relation is valid in wider en- 
ergy bands. In addition, from the results of Hakkila et 
al. (2008) , it is shown that the spectral lag characterizes 
each pulse rather than the entire burst. 

From the theoretical point of view (e.g., Qin et al. 

(2004) ), the rise phase timescale may be responsible for 
the intrinsic pulse width, while the decay phase timescale 
may be determined by geometrical effects (e.g., the curva- 
ture effect). The curvature effect (Qin (2002) , Qin & Lu, 

(2005) , Lu et al. (2006) ) arises from relativistic effects in 
a sphere expanding with a high bulk Lorentz factor V = 
1/(1 - 1 ) 1 ! 2 ~ 100. Because of the curvature of the emit- 
ting shell, there will be a time delay between the photons 
emitted simultaneously in the comoving frame from differ- 
ent points on the surface. However, Zhang et al. (2007) 
showed that the curvature effect alone is not enough to 
explain energy- dependent pulse properties obtained from 
the systematic analysis of lag and temporal evolution. 
Alternative models are the off-axis model proposed by 
Ioka & Nakamura (2001) and the time-evolution of shock 
propagation (Daigne & Mochkovitch (1998) , Daigne & 
Mochkovitch (2003) , Bosnjak et al. (2009) ) may also re- 
produce the spectral lag and the lag-luminosity relation. 
Thus, it is not clear that either the curvature effect or 
other effects cause the spectral lag. While the curvature 
effect should necessarily affect the pulse profile, the time- 
evolution of shock propagation or off-axis model strongly 
depends on unknown model parameters. 

In this paper, in order to unveil the properties of the 
spectral lag for each pulse, we investigate the HETE-2 
sample with a wider energy range especially at the low- 
energy end (>2 keV) than the BATSE sample. In sections 
2 and 3 we explain the sample selection and the pulse- 
fit method. In section 4, we describe the result of the 
obtained relations between the spectral lag and other ob- 
servables, and we discuss a detailed energy dependence for 
the spectral lag in section 5. Finally, we briefly comment 
on the future prospects in section 6. 

2. HETE-2 Sample and Selection 

HETE-2 had two scientific instruments on-board which 
are relevant to our study: the FREnch GAmma-ray 
TElescope (FREGATE), which gave the trigger for GRBs 
and performed spectroscopy over a wide energy range 
(6-400 keV); and the Wide-field X-ray Monitor (WXM), 



which was the key instrument to localize GRBs to ~10', 
and sensitive to the 2—25 keV energy range, lower than the 
FREGATE one. The instruments have two types of data. 
The survey data were recorded with fixed energy bands 
and time resolution whenever the instruments were on. 
The time-tagged data were produced with a fixed duration 
(several minutes) when the instruments were triggered by 
bursts. From the time-tagged data, we can produce light 
curves in arbitrary energy bands, while the BATSE de- 
tector in general created light curves only in fixed energy 
bands (Although the BATSE detector actually has time- 
tagged data, many BATSE GRBs are not fully covered 
due to the limitation of the memory size for the time- 
tagged data.). 

We perform the spectral-lag analysis using a sample of 
8 GRBs detected by HETE-2 with known or estimated 
rcdshifts for the study of the lag-luminosity relation in 
section 4. Our selection criteria for the GRB samples are 
the following: 1) Tgo > 2 s, where Tgo is the observed 
duration including 90% of the total observed counts, and 
2) time-tagged data are available. For the latter, we note 
that the time-tagged data were lost for some bursts due 
to downlink problems or invalidation of the instruments 
(e.g., GRB030328, GRB030329 etc.). For these bursts, 
since the available energy band is too coarse for the sur- 
vey data (e.g., 6-40 keV, 6-80 keV, and 32-400 keV 
for FREGATE), we cannot conduct a detailed study of 
the spectral lag. For the analysis in section 4, we use 
the FREGATE instrument alone because off-axis photon 
events were partially coded and the number of events de- 
tected by the WXM instrument was often small, while the 
FREGATE instrument detected more photons compared 
to those of the WXM instrument due to its relatively large 
effective area (-~150 cm 2 ); not all the selected GRBs have 
enough photons to perform the analysis in the WXM en- 
ergy band. 

In addition, for studying the detailed energy depen- 
dence of the spectral lag for individual GRB pulses in 
section 5, we add 2 GRBs without known rcdshifts having 
sufficiently non-overlapped pulses to the sample. In this 
analysis, we use not only the FREGATE instrument but 
also the WXM instrument, because some GRBs have good 
enough statistics detected by the WXM instrument. Here, 
since there are not good statistics in the multiple energy 
bands for GRB 020124 and GRB 041006, we exclude the 
GRBs from the sample. 

We show the list of 10 GRBs in Table 1 and the energy 
bands in the burst rest frame which are covered by the 
WXM and FREGATE instruments for the selected GRBs 
with known rcdshifts in Fig. 1. 

3. Method 

Each GRB pulse is fitted with a four-parameter pulse 
model (Norris et al. (2005) ), if t > t start 

I{t) = A\exp(-n/(t - i staxt ) - (t - Wt)/T 2 ) + B(K)l) 

= v4Aexp(-Ti/(t + T pcak - tpcak) 

-(i + Tpoak- t pcak )/T 2 )+B(t) (2) 
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Table 1. GRB samples 



GRB 


redshift 


Reference 


010921 


0.45 


Djorgovski ct al. (2001) 


020124 


3.20 


Hjorth ct al. (2003) 


020127 


1.9 1 


Bcrgcr ct al. (2007) 


021211 


1.01 


Vreeswijk ct al. (2003) 


030528 


0.78 


Rau ct al. (2005) 


030725 




Pugliese et al. (2005) 


040924 


0.86 


Wicrsema et al. (2004) 


041006 


0.72 


Stanck et al. (2005) 


050408 


1.24 


Bcrgcr ct al. (2005) 


060121 




de Ugarte Postigo ct al. (2006) 



1 : this is a possible value estimated from the afterglow 
investigation and spectral energy distribution. 

Energy bands in the burst rest frame (WXM & FREGATE) 



> 

S 3 

£ 10' 



te '.<J7, "'0. 



Fig. 1. Energy bands in the burst rest frame for the se- 
lected GRBs. The red dotted bars represent the WXM 
bands and the blue solid ones represent the FREGATE 
bands. The adopted energy ranges are the range be- 
tween the dashed lines (20-100 keV and 100-500 keV). 

and if t < t staY t, I(t) = B(t) , where I is the intensity, t 
is the time after the trigger, r% and t-i are the pulse rise 
and pulse decay constants, A = exp (2(ri/r 2 ) 1/ ' 2 ) , i pca k is 
the time of the pulse's maximum intensity A, i s t ar t is the 
start time, r pca k = (T1T2) 1 / 2 is the peak time from the 
start time i star t, so that t pca k = £ s tart + Tpeak, and B{t) is 
the background function (we utilize a constant or linear 
function). In Eq. 2, i pea k is treated as a primary fitting 
parameter in order to estimate the uncertainty in i pea k 
directly in the fitting procedure. 

The time £ s tart is the formal onset time and in some 
cases t stax t is not indicative of the visually apparent on- 
set time. Especially in the case of t\ ^> 1 s, t stm -t is ex- 
tremely far from the peak of the pulse. Here, as described 
in Norris ct al. (2005) , wc introduce an effective onset 
time i ff arbitrarily defined as the time when the pulse 
reaches 0.01 times the peak intensity. Furthermore, the 
values of t e g are different in different energy bands. For 
HETE-2 GRBs, the statistics of GRBs arc not as good 
as those of BATSE because, e.g., the effective area of the 
FREGATE detector (~ 150 cm 2 ) is lower by a factor of ~ 
10 than that of the BATSE detector (~ 2000 cm 2 ). This 
causes t e g to be scattered in different energy bands due 



to the uncertainties in the determination of t\ and r 2 . To 
avoid this, we adopt an onset time of the "bolomctric" 
light-curve profile, t' cS , derived by fitting the light curve 
in the 6 — 400 keV band, which corresponds to the en- 
tire FREGATE-energy band. The adoption of t' cS is sup- 
ported by Hakkila & Nemiroff (2009) . They showed that 
the onset of GRB pulses occurs simultaneously across all 
energy bands. Thus, we define Xpoak = £ P cak — t' cff in this 
paper. The corresponding uncertainties are calculated us- 
ing the error propagation formula. 

Spectral peak lags are defined as the difference between 
the maximum-intensity times in different energy bands as 



'Hag — ^ P cak,low ^peak,high; 



(3) 



where "low" and "high" represent the low and high en- 
ergy bands, respectively. Another measurable pulse prop- 
erty is the pulse duration w = 3t2(1 



-4v/^7^/3) 1/2 dc- 



fined as the time intervals where intensities are equal to 

^ -^(^ P eak)- 

4. Relation between the Spectral Lag and Other 
Parameters 



In this paper, we adopt two sets of energy bands in 
order to calculate the spectral lag between the two di- 
vided bands. The first one is 6-25 keV and 50-400 
keV in the observer's frame. Although in the previ- 
ous study by BATSE the energy bands between 25—50 
keV and 100—300 keV had been adopted, we adopt the 
lower energy band (<25 keV) and test if the same re- 
lation (e.g., lag-luminosity relation) is established or not. 
Furthermore, for all the previous studies of the spectral lag 
the energy bands refer to the observer's frame. However 
if the spectral lag is a characteristic property of GRBs, 
it is better to derive the spectral lags between the energy 
bands in the burst rest frame. The HETE-2 time-tagged 
data have an advantage for such an analysis, compared 
with the BATSE detector. Thus we adopt energy bands 
20-100 keV and 100-500 keV in the burst rest frame to 
be covered by the FREGATE instrument. The adopted 
energy bands are shown as horizontal lines in Fig. 1. 

Here we use pulses that satisfy the following require- 
ments: the significance of the spectral lag cr s i g > 1.5 and 
positive lag ri ag > 0, where 

CTsig = Ti ag /er lag , (4) 

= (Cpcak.low + CT P o a k,high) • (5) 

The value <7i ag represents the l-cr uncertainty of ri ag and 
c pea k is defined as the l-cr uncertainty of i pea k- For the 
negative lag, its significance is low (eriag < 1) and the num- 
ber of pulses having a negative lag is very small, and we 
do not take it into account. One of the pulse- fitted results 
is shown in Fig. 2 (GRB 050408) using the \ 2 fitting rou- 
tine. For the fitting, we fit pulses to make the obtained 
values of x 2 /d.o.f. (degree of freedom) reasonable (^1). 
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Lag vs L : observer frame 
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Lag vs Duration (50-400keV) : observer frame 
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Lag vs Duration (6-25keV) : observer frame 



to-' 



Lag vs E : observer frame 

peak 



Fio 3 F" 



10H 
10- 2 



10"' 



+ 010921 
▼ 020124 
A 020127 
■ 021211 
• 030528 
□ 040924 
A 041006 
050408 



Lag [si 



* 010921 
T 020124 
A 020127 
■ 021211 

• 030528 
□ 040924 
A 041006 
050408 



Lag [si 



Fig. 3. Scatter plot (observer's frame) of spectral lag vs. luminosity (top left), duration for the low-energy band (top right), duration 
for the high-energy band (bottom left), and -E pea k (bottom right). The dashed lines represent the best-fit functions. 
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Fig. 2. Pulse fit for GRB 050408 in the 6-25 keV and 50-400 
keV bands in the observer's frame. 

4-1. Observer's frame 

First, we show the scatter plots of the spectral lag ri ag 
and luminosity L lso in the top left panel of Fig. 3 where 
the spectral lag is calculated in the observer's frame. The 
luminosity in this paper is defined as the average lumi- 
nosity over the pulse FWHM timescale. As shown in this 
figure, we can see an anti-correlation between the spectral 
lag and the luminosity. Here, for the correlation coeffi- 
cient, we adopt the Spearman rank-order correlation test. 



Furthermore, to estimate the correlation coefficient on the 
basis of the spectral lag's confidence level, we perform a 
Monte Carlo simulation. Since we have already obtained 
Tiag, flag, the luminosity and its uncertainty, we can gen- 
erate a pseudo plot based on the specific probability dis- 
tributions, that is, make a plot similar to the top left panel 
of Fig. 3 with random number seeds. Then we calculate 
the value of the correlation coefficient R for the generated 
pseudo plot using the Spearman rank- order correlation 
test. Finally we repeat the same procedure 10000 times 
with different random number seeds. As we obtain the 
histogram of the correlation coefficient, we regard the 1-er 
width as the 1-cr confidence level. We adopt this method 
in the following analysis. 

For the lag-luminosity relation in the observer's frame, 
we obtain the correlation coefficients as R = -0.79+QQg 
with a chance probability of 7.7 x 10~ 4 at the most prob- 
able value. The best- fit functional form is log(L5i) = 
A x +Bilog(n ag ) with A l = -0.79±0.04,Bi = -1.16±0.07; 
the reduced chi-square is 133.2/12 in the observer's frame. 
Although there are large scatters in the data from the 
best-fit line, the lag-luminosity relation holds even for the 
low energy band (< 25 keV). While our results reconfirm 
the lag-luminosity relation previously reported, our spec- 
tral lag index (-1.2) is slightly smaller than that of Hakkila 
et al. (2008) (index ~-0.6). The slight difference seems 
to come from the following: (1) the different timescales 
to estimate the luminosity (BATSE used 256 ms, while 
we adopt the pulse FWHM timescale), (2) the small num- 
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bers of GRB samples for both Hakkila et al. (2008) and 
HETE-2, and (3) the difference in the adopted energy 
band and/or the instrumental response between Hakkila 
et al. (2008) and ours. Thus, the slight difference in the 
power-law index of the correlations between Hakkila et al. 
(2008) and our results is not surprising. 

We show the scatter plots of the spectral lags and the 
durations at low and high energies in the top right and 
bottom left panels of Fig. 3, respectively. We find corre- 
lations between the spectral lags and durations in both the 
6-25 kcV and 50-400 keV bands. The best-fit functional 
form of these relations is log(w\ ow ) = A 2 + i?2log(ri ag ) with 
A 2 = 1.05 ± 0.06, B 2 = 1.16 ±0.09 and the reduced chi- 
square is 43.0/12 in the 6—25 keV band and log(u>hi g h) = 
A 3 + B 3 log(riag) with A 3 = 0.69 ± 0.05, B 3 = 0.94 ± 0.08, 
and 45.1/12 in the 50-400 keV band. The correlation co- 
efficients are R = 0.66±^ and R = 0.74±°;°^, and the 
corresponding chance probabilities at the most probable 
values are 1.0 x 10~ 2 and 2.5 x 10~ 3 in the 6—25 keV and 
50—400 keV bands, respectively. The low chance proba- 
bilities assure the tight correlations. These results are al- 
most consistent with those of Hakkila et al. (2008) (index 
0.85). 

Finally, the scatter plot of the spectral lag and peak- 
time Speak is shown in the bottom right panel of Fig. 3 
in the observer's frame . The best-fit functional form is 
log(£ P cak) = A 4 + B 4 log(T lag ) with A A = 1.88 ± 0.02, B A = 
-0.31 ±0.02 and the reduced chi-square is 97.8/12. The 
correlation coefficient is R = -0.66io Q8 with chance prob- 
ability 1.0 x 10~ 2 . Thus, we obtain a possible anti- 
correlation between E'pcak and the spectral lag, though 
the dependence of the spectral lag on i^ak is relatively 
weak (index ~ -0.3 ) compared with the other parameters. 

4-2. Burst Rest frame 

We show the results of the relations between the spec- 
tral lag and other parameters in the burst rest frame in 
Fig. 4 (The result of the fitted pulses is shown in Fig. 
8), and the adopted energy bands are determined to have 
the same energies in common (20—100 keV and 100—500 
keV). As is the case for the observer's frame, the best- 
fit parameters of each relation are summarized in Table 2. 
We find that there is no significant difference in results be- 
tween the observer's frame and the burst rest frame. The 
obtained results support the idea that the cosmological ef- 
fects should not significantly change our measurement in 
the observer's frame even though most GRBs are found 
at high redshifts. The intrinsic properties for GRB pulses 
predominate over the cosmological effects as suggested by 
Hakkila et al. (2008) and Hakkila & Cumbee (2009) . 

4-3. Discussion 

We have obtained the correlation between the spectral 
lag and duration, L lso and £p ak in the observer's and the 
burst rest frames. In particular, our result extends the 
energy coverage to a lower energy band (6—25 keV). This 
indicates that the GRB emission in the wide X-ray band 
has the same origin. 

As there is no significant difference between the results 



in the observer's and burst rest frames, it is natural to 
adopt the burst rest frame to discuss the origin of the 
spectral lag. Thus, in the following discussion, we refer to 
the case of the burst rest frame. 
4-3.1. Physical Origin of the Relations 

To account for the relation between the spectral lag and 
Liso, let us consider the off-axis model suggested by Ioka & 
Nakamura (2001) ; the detector observes a GRB jet with 
different viewing angles 9 V . The intrinsic physical param- 
eters, bulk Lorcntz factor T, opening half-angle A9 ~ 1/r, 
shell radius ro from the center, and £p oak in the comoving 
frame are assumed to be the same for all GRBs. In this 
paper, we adopt the same parameters as those of Ioka & 
Nakamura (2001) ; TA9 = 1, r Q /c/3T 2 = 1. This model 
assumes an intrinsic spectral shape in the comoving frame 
which is approximated by the Band function (Band et al. 
(1993) ) as 



f(E') 



£1 



l+a B 



1 



f>B-"B 



(0) 



where as and /3b are the low- and high-energy indices, 
s describes the smoothness of the transition between the 
high and low energies and Eq is the break energy. Ioka 
& Nakamura (2001) showed that the observable values 
such as luminosity, the pulse duration (FWHM) w, and 
the peak energy -E pca k7 change with the viewing angle 9 V 
and correlate with the spectral lag as, 



Lisa oc r lag s+1 , 

1 

w oc 1 + const x r^g 1 . 



peak 1 



oc 1 + const x r, 



(7) 
(8) 
(9) 



We have obtained L; so oc r lag ' from the lag-luminosity 
relation in Fig. 4 so that the off-axis model with s = 
1.5 and the typical value for the low-energy photon index 
as = — 1 can reproduce the lag-luminosity relation. 

The expected theoretical results are superimposed on 
Fig. 4. Here, we adopt arbitrary normalization values for 
the theoretical lines. For the lag-duration relation in the 
high-energy band (100—500 kcV), the observational points 
are consistent with the theoretical curve. For the lag- 
duration relation in the low-energy band (20—100 keV), 
the observational points agree with the theoretical curve 
for small spectral lags, although there are some outliers for 
large spectral lags. For the lag-E? pea k relation, we find a 
consistency between the observational points and the the- 
oretical curve. Except for some outliers, we find that the 
off-axis model can explain the observational results well, 
even though the model seems to be oversimplified. (All 
intrinsic physical parameters are common for all GRBs 
in this model.) Furthermore, for the off-axis model the 
spectral lag is calculated using the difference between the 
peak times in the different energy bands just as we calcu- 
late the spectral lag, unlike Norris et al. (2000) , in which 
the spectral lags are calculated using the CCF method. 
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Lag vs Duration (20-100keV) : rest frame 
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Fig. 4. Scatter plots (burst rest frame) of spectral lag vs. luminosity (top left), duration for the low-energy band (top right), 
duration for the high-energy band (bottom left), and E pea k (bottom right). The dashed lines represent the best-fit functions. For 
the last three plots, the long dashed lines show the theoretical curve expected in the off-axis model (Ioka & Nakamura (2001) ). The 
dotted lines represent the best-fit functions including the systematic uncertainty and excluding GRB 030528. 



4-3.2. Yonetoku Relation 

We now consider the consistency with the Yonetoku 
relation (Yonetoku et al. (2004) ) in this section. In 
the preceding section, we have found that the off-axis 
model reproduces the observational results and have not 
imposed any limitations such as the Yonetoku relation 
(L iso oc£> ak ). 

Assuming that the Yonetoku relation is valid, from our 
result on the lag-luminosity relation ( ii SO oc Tj~g' 23±0 ' 07 ), 
the lag--Bp Ca k relation is expected to satisfy 

r l/2 -1.23/2 -0.62 / 1n \ 
-Epeak,cxp OC L is ' CXT lag OC Tj ag (10) 

The index (-0.62) is small compared with the obtained 
result (-E P cak,obs oc t,""' 33 * - 03 ). Note that the subscripts 
"exp" and "obs" represent the expected and observed val- 
ues, respectively. 

Let us assume that the determination of the spectral 
lag has a systematic uncertainty <7 sys of 0.05 s resulting 
from the overlaps of the GRB pulses or some calibration 
uncertainties. In addition we exclude the peculiar case 
of GRB 030528, which has a very long spectral lag that 
may be due to the overlaps of multiple pulses. Then, 
the best-fit functions become log(L.5i) = (—1.38 ±0.13) — 
(1.59±0.28)log(7i ag ) with reduced chi-square xl = 8.2/12, 
and log(-B pcak ) = (1.63 ±0.08) - (0.76 ± 0.17) log(r lag ) with 
reduced chi-square xl = 12.5/12. The Yonetoku relation 
and the revised lag-luminosity relation give us 



-C'peak.exp OC T lag oc T lag ) K 11 ) 

which agrees with the revised result (index —0.76). 
Therefore, considering the small sample and observa- 
tional uncertainties, we cannot exclude the validity of the 
Yonetoku relation from our results. 

5. Detailed Energy Dependence of Spectral Lag 
and Other Properties 

Next we consider the detailed energy dependence of the 
spectral lag and other properties (the durations includ- 
ing the rise and decay times), besides the lag- luminosity 
relation described in the former sections. 

5.1. Energy Dependence of Duration, Rise and Decay 
Phase 

Zhang et al. (2007) studied the energy dependence of 
temporal properties represented by the formulae, 

woc£ Q ", (12) 
T lisc = ^{w-r 2 )KE a "°°, (13) 

Tdecay^(u> + T 2 )cx£; ad -, (14) 

where r r i S0 and Tdecay are the rise and decay timescales 
(Norris et al. (2005) ). They found that a w and ayee 
are highly correlated, while a w and a r i se are not strongly 
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correlated. Here, it may be reasonable to assume that 
the intrinsic pulse width is responsible for the rise phase 
timescale, while the decay phase timescale is determined 
by the geometrical effect in the relativistic expanding 
shell. Furthermore, the decay time interval dominates the 
duration because the typical pulse shape shows a fast rise 
and exponential decay (FRED). Thus, it is natural that 
the decay phase is highly dependent on the duration and 
the rise phase is not strongly related to the duration (or 
decay time). 

We show the results of the plots among a w , a r j S e and 
adcc in Fig. 5 (The result of the fitted pulses is shown 
in Fig. 9). The top panel of Fig. 5 shows the scatter 
plot of a w versus a r ; sc in our HETE-2 sample. Although 
the uncertainty for each point is very large, we find a 
marginal linear relation with correlation coefficient R = 
0.51±g;||. The best-fit function is a risc = (0.04 ±0.15) + 
( 1.06 ± 0.41 )a w . If the rise timescale is determined only by 
the intrinsic pulse width, the correlation between a r i SG and 
a w should be weak. However, from this result a w seems 
to be roughly proportional to a r i S0 . Therefore, the rise 
timescale depends not only on the the intrinsic pulse width 
but also somewhat on the geometrical (curvature) effect. 
Some previous studies may give us clues to understanding 
this behavior; Lu et al. (2007) and Peng ct al. (2009) 
found that -E p0 ak decays monotonically through long GRB 
pulses. This energy decay occurs even prior to the pulse 
peak, namely in the rise time phase. Therefore, the pulse 
rise phase is a part of the .©peak-decay phase. Hakkila & 
Cumbee (2009) also demonstrated that the high-energy 
pulse intensity is starting to decline prior to the pulse 
peak in the low-energy band, as is the case for the HETE- 
2 results. Thus, these results indirectly indicate that the 
pulse rise timescale is affected by the pulse decay time. 

The middle panel of Fig. 5 shows the scatter plot of 
a w versus ckdec- The best-fit function is adcc = (—0.01 ± 
0.13) + (1.03 ± 0.37)a w with correlation coefficient R = 
0.671q;37. Since a relatively good proportionality between 
a w and ayec exists, this result supports the approximation 
w ~ r dccay and the assumption that the curvature effect 
determines the decay timescale. 

The relations between a r i so versus adcc are plotted in 
the bottom panel of Fig. 5, where the best- fit function is 
a do c = (-0.04 ± 0.13) + (0.99 ± 0.40)a risc with correlation 
coefficient R = 0.46lo'3g. The result also seems to show 
that the rise timescale is slightly affected by the curva- 
ture effect. Although the uncertainties in the correlation 
coefficient and fitting parameters arc large, the relation- 
ships between a w , a r i se and adec are consistent with those 
of Zhang et al. (2007) (the three functional forms agree 
with ours). 

Shen et al. (2005) computed the temporal profiles of 
the GRB pulse in the four BATSE energy bands, with the 
relativistic curvature effect of an expanding shell. They 
included an intrinsic "Band" shape spectrum and an in- 
trinsic energy-independent emission profile, and estimated 
the dependence of the duration and other properties on en- 
ergy as w cx iT - 2 — - 1 (a w ~ a dcc = -0.2 to -0.1). On 
the other hand, Daigne & Mochkovitch (1998) calculated 



the time-evolution of the internal shocks (hydrodynami- 
cal effect) assuming a highly non-uniform distribution of 
the Lorentz factor, and obtained the energy dependence 
as w oc E~ 0A . 

In our result, shown in Fig. 5, a w and adcc range from 
-0.8 to and the expected energy dependences for the 
models of Daigne & Mochkovitch (1998) and Shen et al. 
(2005) are represented as a long dashed line and a shaded 
portion, respectively. The data points are widely scattered 
so that the simple model of Daigne & Mochkovitch (1998) 
or Shen et al. (2005) alone cannot explain the results we 
obtained. 

5.2. Physical Origin of the Spectral Lag for Individual 
Pulses 

In this section, we try to clarify the origin of the spectral 
lag of our HETE-2 GRBs, apart from the lag-luminosity 
relation described in the preceding section. First we need 
to check whether only the simple curvature effect, which 
should be in any case included, can explain the energy 
dependence of the lag or not. Here we consider the cur- 
vature effect model described by Lu et al. (2006) . They 
calculated the light curves from an isotropically expanding 
sphere with a constant bulk Lorentz factor and the Band 
function for a rest-frame radiation spectrum. A Gaussian 
pulse was assumed for the light curve in the source rest 
frame. While the generic formula for the spectral lag is 
complicated, they demonstrated that the spectral lag has 
a energy dependence with lag oc E below a saturated en- 
ergy, E s = 1.67£'pcak for the typical parameter sets of the 
low-energy index cib = -1, high-energy index (3b = -2.25 
and the shell radius ro = 3 x 10 15 cm. Considering the 
beaming effect, photons of energy at E > E s would mainly 
come from the area of the surface around the line of sight, 
i.e., < r _1 (where 9 is the angle to the line of sight). 
When E > E s , the contribution to the corresponding light 
curve largely comes from the high-energy portion of the 
rest frame spectrum, which causes the peak time of the 
light curve to change less, and the lag would saturate. On 
the other hand, for E < E s , "off-axis" (6 > T^ 1 ) photons 
may contribute to the light curve so that the lag increases 
with the increasing energy difference in two energy bands. 

We choose the peak time at the lowest energy (~ 1 
keV) arbitrarily to match the observational value of T pea k. 
Lu et al. (2006) showed that the lag does not depend 
strongly on the radius ro- So based on their results, we 
write approximately 



Tla 



where a = 10 



aE, 
aE s 

r 

Too 



if E < E s 
if E s < E 

-2.8 

[s/keV]. 



Then Tl 



peak 



is written as 



T, 



peak 



-aE, 
- aE. 



itE< E s 
\fE,<E 



(15) 
(16) 

(17) 



where t is the peak time at the lowest energy. Using Eq. 
17, we try to reproduce the spectral lag for the examined 
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Fig. 5. The scatter plots for indices of o w , a r ise and a^ ec . 
(Top): a w vs. a r i so . (Middle): ct w vs. ayec- (Bottom): 
a riso vs - Q dcc- The dashed line shows the best-fit lin- 
ear function. The shaded area and long vertical and hori- 
zontal dashed lines represent the expected values from the 
simple curvature (Shcn et al. (2005) ) and hydrodynam- 
ical effects (Daigne &; Mochkovitch (1998) ), respectively. 



pulses with the curvature effect by adjusting bulk Lorentz 
factor r in Fig. 6 and 7. Here, the energy and T poa k 
are translated into the burst rest frame ones with known 
redshifts, and for GRB 030725 and GRB 060121 with- 
out known redshifts assuming that their redshifts are 1. 
Although our fits are based on an empirical formula with 
a particular parameter set (olb = -1, Ab = -2.25), the data 
points do not largely contradict the tendency predicted by 
the curvature effect at a particular bulk Lorentz factor T, 
as shown in Fig. 6. But for some pulses such as GRB 
021211, the 2nd pulse of GRB 040924, the 2nd pulse of 
GRB 050408 and GRB 030725, the model we adopt cannot 
reproduce the spectral lag well as shown in Fig. 7. Even 
for such pulses changing the parameters that we have fixed 
here may soften the contradictions. Alternatively, the off- 
axis model or the temporal evolution of the internal shock 
propagation may play an important role in the spectral 
lag, or a pulse-overlap effect may be included. 

To examine whether the pulse duration and spectral 
lag arc explained synthetically by the curvature effect, wc 
summarize the results of the estimated T, a w and other 
properties in Table 3. The empirical formula based on Lu 
et al. (2006) has been derived from the assumption for 
the intrinsic pulse duration At mt = 10 5 s in the source 
frame. Although there is no reason to adopt this value, 
the resultant timescales estimated from the obtained T 
are of the same order of magnitude as 10 5 s. Even for the 
GRBs whose spectral lag can be explained by the curva- 
ture effect, it is hard to confirm the consistency between 
the prediction of the energy dependence of the pulse du- 
ration by the curvature effect (a w = -0.2~-0.1) and the 
experimental values because of the large uncertainties in 
a w . However, we may say that the model based on the 
curvature effect does not contradict both the spectral lag 
and duration at a particular bulk Lorentz factor T. Thus, 
the spectral lag can be a tool to help estimate the bulk 
Lorentz factors. 

Since in this analysis only a finite energy range (2—400 
kcV) is available, there are only a small number of points 
for the energy range, where the lag is saturated above E s . 
For many GRBs, wc could only plot one point above E s , 
which leaves the possibility that a significant lag takes 
place above E s (no saturated energy). Even in the study 
by Liang et al. (2006) , although the peak energy Speak 
is ~ 54 keV, they also could plot only one point for T pca k 
above Speak due to the poor effective area for the higher 
energy ranges. To clarify the origin of the spectral lag 
further, wc need to detect GRB photons in the higher 
energy ranges above S pea k to describe the light curve and 
determine T p0 ak with confidence. 

6. Towards a Unified Theory 

Our results suggest that there are correlations between 
the spectral lag and other observational properties for 
each GRB pulse. The correlations we found and some 
additional pulse correlations such as spectral hardness or 
pulse asymmetry etc. reported by Hakkila & Cumbce 
(2009) are important hints to specify or constrain models 
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Fig. 6. Energy vs. ^peak plots in the burst rest frame with the theoretical model lines, showing good agreement 
with the curvature case. Each line represents the curvature-effect line at the corresponding bulk Lorentz factor T. 
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Fig. 7. Energy vs. T poa ]< plots in the burst rest frame with the theoretical model lines, showing poor agreement 
with the curvature case. Each line represents the curvature-effect line at the corresponding bulk Lorcntz factor T. 



Table 2. Fitting parameters of the lag-luminosity -duration -i? pGa k relations 







A 


B 


reduced chi-square 


correlation coefficient 


luminosity 1 


obs 


-0.79±0.04 


-1.16±0.07 


133.2/12 


-0 7Q+ U i6 

u -' a -0.05 




rest 


-1.09±0.04 


-1.23±0.07 


97.1/13 


-0.90iS 


duration 


obs 


1.05±0.06 


1.16±0.09 


43.0/12 




(low energy) 


rest 


1.06±0.04 


1.15±0.06 


72.5/13 


n 74+ 06 

u - '^-0.15 


duration 


obs 


0.69±0.05 


0.94±0.08 


45.1/12 


u - '^-0.14 


(high energy) 


rest 


0.67±0.06 


0.73±0.09 


26.6/13 


72+ 07 

u - ' z -0.22 




obs 


1.88±0.02 


-0.31±0.02 


97.8/12 


-o.66±8:S 




rest 


1.82±0.02 


-0.33±0.03 


33.9/13 


-o-si±8-iS 



; — u.uo 

: log(L 5 i) = A + Blog(n ag ), °: log(w) = A + Slog(n ag ), *: log(£:p 0ak ) = A + _Blog(Ti ag ) 



Note that "obs" represents the observer's frame and "rest" represents the rest frame. 



of GRB prompt emission. 

On the other hand, because of the large dispersions, 
the spectral lag relations are not so useful as tools to 
measure cosmological distances so far, compared with the 
Yonetoku relation. We should note that there may still 
be systematic uncertainties in the obtained lags, which 
may change the correlations as discussed in §4.3.2. While 
the obtained lag-luminosity, --E poa k and -duration rela- 
tions can be consistent with a specific model, namely the 
off-axis model suggested by Ioka & Nakamura (2001) , the 
energy dependences of the spectral lag seem to be con- 
sistent with the simple curvature effects for some GRB 
pulses. The assumptions inferred in Ioka & Nakamura 



(2001) and Lu et al. (2006) arc different so that we have 
discussed the correlations and energy-dependences in the 
spectral lag with two independent models. Although such 
methods do not give us a consistent picture for the spec- 
tral lag so far, the discussion in this paper may help to 
determine which models are more appropriate. 

For a unified theory to explain the spectral lag and other 
temporal spectral characteristics, the effect of the curva- 
ture, viewing with an offset angle to the jet, time-evolution 
of shock propagation, and other effects must be taken into 
account synthetically and theoretical investigations need 
to be done. To have further quantitative discussions, we 
need a sample which includes many GRBs having a good 
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Table 3. Summary of the obtained bulk Lorcntz factor, the power-law index of the duration, the observed duration and the expected 
intrinsic pulse width 



GRB 


r 




7/7 U [sl 


Ar xp s 


010921 


60 




21.4±2.9 


0.77x10 s 


020127 


500 


-0 67±0 34 


70+0 01 


1.75xl0 5 


030528 (1st pulse) 


50 


-0.63±0.61 


26.1±14.7 


0.65xl0 5 


030528 (2nd pulse) 


30 


-0.42±0.15 


58.8±6.8 


0.53xl0 5 


040924 (1st pulse) 


140 


-0.18±0.47 


1.8±0.5 


0.35xl0 5 


050408 (1st pulse) 


120 


-0.20±0.30 


2.8±0.2 


0.40xl0 5 


060121 


160 


-0.42±0.17 


1.8±0.2 


0.46 xlO 5 


021211 




-0.58±0.14 


2.4±0.1 




030725 




-0.18±0.07 


19.7±0.5 




040924 (2nd pulse) 




-0.51±0.19 


1.4±0.1 




050408 (2nd pulse) 




-0.22±0.36 


2.0±0.2 





Note that the bulk Lorentz factor T is estimated by the lag analysis in Fig. 6 and w bs is the observed duration in the 
burst rest frame and Ar^* t p is the expected intrinsic pulse width from the obtained bulk Lorentz factor F. 



S/N ratio detected in a wide band (keV— GeV) with ob- 
servationally known redshifts. 
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Fig. 8. The result of the fitted pulses in the burst rest frame (20-100 keV and 100-500 keV). 
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Fig. 9. The result of fitted pulses in multiple energy bands. 



